Energy and information flow in superlensing 

E. Fourkaj*) and I. Velchev 
Department of Radiation Physics, Fox Chase 
Cancer Center, Philadelphia, PA 19111, U.S.A. 

A. Smolyakov 
Department of Physics and Engineering Physics, 
University of Saskatchewan, Saskatoon, Canada 
(Dated: March 4, 2009) 

Abstract 

In the superlens problem, the presence of the reflected wave plays a critical role in the energy 
and information transfer from the object to recording device: both incident and reflected evanes- 
cent waves are required to insure nonzero energy /information flux from the object to the recorder. 
Therefore an optimization is required between the image quality (characterized by a transfer func- 
tion) and information transfer rate (characterized by the energy/information flux). It is shown 
that the introduction of the recording device may lead to the deterioration in the image quality as 
compared to the ideal case when no recording device is present. The decline in the image quality is 
due to the generation of phase-shifted reflected evanescent wave between the lens and the recorder. 
It is also shown that the magnitude of the energy flux depends on the amount of dissipation in the 
lens and the recorder in a non-monotonous way. 
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I. INTRODUCTION 



Propagation of the electromagnetic radiation in materials with negative dielectric per- 
mittivity (and permeability) or the so called left-handed materials (LHM) has attracted 
great deal of attention in recent years [l|, |2j, y|, |4|. The increased interest in the properties of 
such media has been driven by their potential applications in various branches of science and 
technology, ranging from development of terahertz computer chips and next generation of in- 
tegrated opto-electronic devices to ultra-sensitive molecular detectors (used in bio-molecular 
imaging applications), and invisibility cloak j^, [(], Q]. One additional application of these ma- 
terials which has been widely studied in recent years is related to the possibility of creating 
the so called superlens: a sub-wavelength optical imaging system beyond the diffraction 
limit 8|. 

The superlens phenomenon is essentially based on amplification of evanescent compo- 
nents of the object spectrum (these are the components that carry sub-wavelength infor- 
mation about an object), facilitated by the excitation of surface plasmons si, 0]. This 
part of the signal spectrum is normally lost in the standard optical devices resulting in 
the diffraction limit. A simplified one dimensional object can be represented in the form 
F(y) = J Fkexp (ik y y) dk y . The components with k y > u/c constitute the evanescent part 
of the spectrum and the components with k y < u/c are propagating. In a standard optical 
device, only the propagating components carry information from the object plane to the 
image plane. The superlens on the other hand, has a property that both propagating and 
evanescent parts of the spectrum can be transmitted via such a lens. The transmission of 
the evanescent part of the spectrum becomes possible via the amplification of the evanescent 
waves due to the resonant excitation of the surface modes. 

The transfer of information costs energy lOUllj. As shown by Bekenstein 10[], there exists 
a universal bound on the energy cost per bit of information transferred E/I > h\n[2]/7iT, 
which means that larger information transfer rates ultimately require higher energy transmis- 
sion (it also confirms the fact that it is impossible to transfer information without spending 
some energy [12]). Therefore, this necessarily means that evanescent waves carrying the in- 
formation about sub-wavelength structure of an object have to transfer energy. Depending 
on the "complexity" of an object (percentage that evanescent components take in its Fourier 
spectrum) , different amount of energy flux is needed in order to transmit this sub- wavelength 
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information from an object plane to an image plane. There is an important difference be- 
tween the energy transport by the propagating (k y <u/c) and evanescent (k y > uj/c)) parts 
of the spectrum. For propagating waves, the energy flux is finite and can be calculated 
independently of the receiver. In other words, a finite energy flux can be set up in a half 
infinite region so that an ideal recording device can be imagined as that having a black body 
type absorption properties, receiving all of the available energy/information. In the evanes- 
cent part of the spectrum however, the situation is radically different. The energy flux in 
a half infinite region (where there is only one evanescent component exp (—ax) present) is 
identically zero. The second component exp (ax) with a proper phase shift is required to set 
a finite energy transport [3] . The existence of the second component and the value of the 
phase shift are intrinsically related to the properties of the recording device. For example, 
in the absence of absorption in the recording device the energy flux in the evanescent part 
of the spectrum is identically zero. Thus, the recording device becomes an integral part of 
the imaging system in the process of capturing the energy and information, carried by the 
evanescent waves. 

Typically the quality of the lens is characterized by its transfer function (an ideal imaging 
system would have the transfer function equal to unity for all k y , see below in the Section 2). 
The transfer function alone however ignores the issue of the transmitted energy flux (hence 
the amount of transmitted information). This is not important for propagating waves as 
the closeness of the transfer function to unity is the only criterion that determines the 
quality of the lens. However, for the evanescent part of the spectrum, the ideal transfer 
function(r = 1) does not guarantee sufficient energy /information flux. Therefore, for the 
subwavelength imaging the optimization of the energy flux (determined by the dissipation 
in the recording device) becomes necessary. Discussion of these issues is a subject of the 
present paper. 

A generic imaging system consists of three main components: (1) the source of infor- 
mation carrying medium (electromagnetic waves, sound waves, etc.), (2) the lens, whose 
role is to apply proper phase adjustments to different spatial components of the spectrum 
and (3) the recording device, whose role is to capture the information, carried by mediating 
waves. By its purpose, the recording device is necessarily dissipative since it has to imprint 
the arriving information into "itself. The presence of the recording device also creates a 
reflected wave that ultimately ensures the energy and information flow between the lens and 
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rd. 

In many previous research studies 1 1 51 ] concerning image formation by the so called 
Pendry lens (a lens made out of a material having negative permittivity but positive perme- 
ability) the authors investigated lens's optical transfer function as well as its ability to focus 
an image in the absence of the recording device. This case however describes unrealistic 
situation, since therein obtained optical transfer function does not contain the influence of 



the recording device. As shown in Ref. 



161 ]. the introduction of the recording device may 



in some cases significantly deteriorate the quality of a resultant image (see below), which is 
due to the generation of a reflected wave by the lossy recorder. 

Here we investigate the role of the IRD (Image Recording Device) in the energy and 
information transfer between the object and image planes. We show that the level of image 
deterioration depends on an interplay between the amount of information needed to describe 
the imaged object (in a given time period) on one hand and the optical transfer function of 
the lens on the other. In particular, we study the case of a silver film (negative permittivity) 
deposited on another medium (with positive permittivity) behind which a recording device 
is present. This particular configuration is descriptive of the near- field optical lithography 
imaging system used in recent investigations [h 17] 



II. OPTICAL TRANSFER FUNCTION OF THE SYSTEM 

The imaging problem can be described in terms of the optical transfer function t(x, k y , u) 
{k y designates the in-plane wave vector directed along the surface of the material), defined 
as the ratio of Fourier components of image field to object field, B^ g (x) / -B o y(0) (for — oo < 
k y < oo) at a given imaging plane x. The transfer function r can be used to find the 
reconstructed field in the image plane in the form 

B img (x,y,t) = J B (ky)r(x,ky,uj)e t{k y y ~^dky, (1) 

where B (k y ) is the wave vector spectrum of the source (imaged object). Thus, the ability 
of the system to image the object is completely determined by the optical transfer function, 
which in itself depends on many physical parameters of the system. In an ideal case, the 
transfer function should transfer all spatial harmonics equally or r(x, k y , k ) = const, for 
— oo < k y < oo. In reality however, the transfer function is a non-monotonous function 
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of the wave vector k y , medium material type, its thickness d, and the position x of the 
imaging plane relative to the position of the object plane. It can be found by considering a 
p-polarized wave (electric vector in the plane of incidence) incident on a thin slab of thickness 
d, dielectric permittivity e\ < and magnetic permeability \i\ = 1 as shown in Figured) The 
lens is separated from the source plane by the medium with permittivity eo and permeability 
fio = 1 and the thickness a. The detector is modeled as a medium with permittivity €3 and 
permeability fa = 1. It is located at a distance c from the lens. The space between the 
detector and the lens is filled with medium having permittivity e 2 and permeability fa — 1. 
The optical properties of this imaging system are obtained by taking the ratio of the field 
in the region x > a + d + c to that at the object plane (in current calculations the object 
plane is assumed to be at x = 0). The electromagnetic fields in each region of interest are 
found from solving the well known wave equation, 

d fldB z \ u 2 / tfc 2 ' 
dx \ e dx ) c 2 V uj 2 



^(-^)+^U--V)^ = 0> ( 2 ) 



with a general solution having the following form, 

B z = (A x e ikx + A 2 e~ ikx ) j^v-**) (3) 



where k = 00 / c^J (e — k 2 c 2 / uj 2 ). The electromagnetic fields in the first medium [x < a) 
represent a sum of incident (emitted by the object at x — 0) and reflected waves. The field 
in the detector (x > a + d + c) is modeled as a transmitted wave only. Matching solutions 
at different boundaries by requiring the continuity of B z and l/edB z /dx across interfaces, 
one arrives at the expression for the transfer function at x = a + d + c (imaging plane), 

g e i(k3C+k!d+koa) t£t 

r(a + d + c, k v , uj) = — . . — (4) 

where E = (£ - - 6) (6 + 6), S x = (Co + 6)(6 - + 6), 2 2 = (Co - 

- e 3 )(6 + £3), S 3 = (£ + £i)(£i + £ 3 )(6 + £3) and £ = A*/^. The zeros of the 
denominator in Eq. H] define dispersion relation for the eigenmodes supported by the given 
system, which in general are coupled plasma surface waves running on either side of the 
lens. Their presence is detrimental for imaging purposes, since they disproportionately 
enhance the resonant spatial frequencies in the image. Presence of dissipation in the lens may 
dramatically reduce and widen the resonances, leading in certain cases to some improvement 



in the image quality 



, [lsj]. As shown in Ref. [3], a favorable distribution of material in 



the asymmetric superlens system should be such that Ree2 = Ree^ = —Ree\. This choice 
of permittivity distribution leads to weaker dependence of the transfer function on the wave 
vector as well as reduces the amount of the reflected wave generated between an object and 
the lens (the presence of large reflected wave leads to distortion of the object field, thus 
producing artifacts in the image). Figure ([2]) shows the reconstructed image of two slits 
of width 2.5 ran each, separated by a distance of 20 nm, for two different cases, with and 
without the recording device. It should be noted here that the absence of the recording 
device corresponds to the case when lme 3 =0 (Re[e 2 ]=Re[e 3 ]). As one can see, the presence 
of a dissipative recording device somewhat deteriorates the quality of an image as compared 
to the ideal case of no detector present. It is also worth noting that the quality of the 
resultant image is not only dependent on the optical properties of the imaging system but 
is also linked to the imaging object itself as well as its spatial spectral characteristics or 
information content. As mentioned earlier, the ideal imaging system should transmit all 
spatial components equally. In reality however, for all practically accessible lens parameters 
(its width, amount of dissipation, and the choice of material surrounding the lens), the 
optical transfer function is usually a non-monotonous function of the wave vector k y with 
zero asymptotic value at large k y . Therefore, the quality of the reconstructed image depends 
on the relation between the characteristic frequency kfjj- of the imaged object (maximum 
value of the wave vector k y for which the object's Fourier spectrum is not negligible) and 
the characteristic frequency k^f of the lens (maximum value of the wave vector k y for which 
the optical transfer function is not negligible) as well as the functional shape of the optical 
transfer function. The condition k°^ < k^jj- along with the optical transfer function having 
weak dependence on k y presents a favorable case in terms of the quality of the resultant 
image. As mentioned earlier, the presence of the recording device reduces the value of 



161 ]. but if the above condition is still satisfied, the image deterioration due to the 
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K eff 

presence of the recording device will be minimal. 



III. INFORMATION CONTENT OF THE IMAGED OBJECT 



In order to quantify the " complexity^ of an object or its information content, we'll use the 



ideas from the communication theory 121]. specifically the definition of information through 
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entropy, 

i 

where are the a priori probabilities of occurrence of various states related in our case to 
the imaging object. It is natural to define it as the probability of occurrence of k l y spatial 
frequency in the Fourier spectrum of the object. In this respect, the object can be viewed 
as an information source which contains a message in a form of a sequence of different 
components of the wave vector k y . This probability can be expressed in a discrete form as, 

rkl+S/2 

p i =p(k t y ,5)= / P(k y )dk y (6) 

Jk*-S/2 

where 5 is the sampling size and P is the density probability function, which is related to 
object's spectral density distribution, 

poo 2 

P(ky) = A 



B,{y)e ik ^dy 



(7) 



where A is a normalization factor found from the condition P(k y )dk y = 1. Substituting 
Eq. [7] into Eq. [6] and subsequently Eq. [5j one obtains an expression for information measure 
1(5) of the imaged object (in bits) for the given degree of discretization 5. Taking subse- 
quently the limit 5 — > provides the total information contained in the imaged object (the 
existence of such a limit follows from the normalization condition on the density probability 
function). If one is able to transfer all of this information from the object plane to the image 
plane and imprint it into the recording device, a perfect image will be formed. However, 
since the recording device in its design is a discrete system (degree of discretization is deter- 
mined by the pixel size), the recorded image will never contain all of the initial information 
"emanating" from the object even in the presence of an ideal lens. Figure [3] shows the de- 
pendence of the information contained in the object (two slits of width 50 nm separated by 
the distance 100 nm) as a function of the sampling size S (the wave vector's sampling size 
5 is related to the recording device's pixel size a through the following relation 5 = 27Tp-, 
where A is the characteristic wavelength of the object's spatial inhomogeneity; it is the dis- 
tance between the two slits in the example presented above). As 5 increases, the information 
imprinted in the device decreases, reducing the image quality of the original object. This 
information reduction effect is due to a coarse-graining or "zooming out" procedure irrevo- 
cably introduced by observing the "world" at a finite level of resolution. The act of recording 



will raise the entropy of the recording device by at least the amount equal to that of the 
recorded information. The entropy of the observed object however should decrease by the 
same amount as a result of the act of observation 19J. This entropy reduction of an object 



can be also understood from considering the relation between the entropy and information: 
the thermodynamic entropy is an estimate of the amount of further information needed to 
define the detailed microscopic state of the system. Once the observer has recorded a part 
of the total information content of an object, the remaining not yet observed or hidden 
information determines the object's entropy 20]. 

In the arguments presented above, we didn't account for the influence of the superlens 
on the object information. Since the lens in general, does not equally transfer all spatial 
harmonics k y but rather enhances one part of the original spectrum and suppresses the 
other (except for an ideal case of extremely thin lens of thickness less than 1 nm, in which 
case the transfer function is indeed constant over very large wave numbers), the object's 
spectral distribution "emanating" behind the lens is somewhat reshuffled (changing the 
probabilities of occurrence of k y spatial components), distorting the shape of the original 
object and changing its information content (Eq. [5]). The object's information reduction- 
distortion effect by the lens on one hand is due to already mentioned presence of surface wave 
resonances. On the other hand it is also due to the so-called "information hiding" effect, in 
which a part of the original object's information imprinted in the evanescent waves becomes 
hidden to the outside world as an "observer" moves away from the object's plane. This 
information loss is reversible in its nature, since the hidden information can be recovered by 
simply putting the recording device closer to the object. The presence of the lens simply 
allows recovering of some of the hidden information at a distance from the object (not 
possible in the absence of the lens) . 

In addition, when there is dissipation in a lens, part of the object's information gets irre- 
versibly lost to the observer due to energy absorption. This information does not disappear 
to the outside world however, but is carried away by the thermal photons radiated away 
from the lens. 
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IV. ENERGY AND INFORMATION TRANSPORT IN SUPERLENSING 



It would be instructive to first revisit the question of what the act of information recording 
is. During this process, the incoming information is imprinted in the recording device via 
some process of energy deposition. Thus the recording device registers the energy flux, or 
the normal component of the Poynting vector but not the value of the field. Therefore, the 
optical transfer function that has been calculated as a ratio of the field value at the imaging 
plane to that at the object plane (Eq. H]) and explored in many previous investigations serve 
more of an instructive role on the nature of the evanescent field amplification rather than 
describing the real experimental situation. Therefore, the experimentally relevant optical 
transfer function should be defined as that proportional to the normal component of the 
Poynting vector at the imaging plane. Using the equations from the previous section, the 
normal component of the time averaged Poynting vector P x = ^Re[E y B* z ] is given by the 
following expression, 

P x = l -Re[ f J ^(k' y )r(k' y )B (k' y )T*(k' y - k y )B*(k' y - k y )e- ik ^dk' y dk y ] (8) 

For the case when the incident wave is a plane wave, so that B (k y ) = B 5(k y — k y ), the 
x component of the Poynting vector simplifies to P x = -Re^s] |r| 2 |-Bo| 2 - Thus, the optical 
transfer function in this case may be defined as, 

T(a + d + c,k° y ,u J )=[\r\ 2 Re[^]} k=k0 . (9) 

y y 

It is interesting to note that the energy flux optical transfer function T is equal to zero when 
-Re [£3] = 0, which for the evanescent part of the spectrum occurs when Jm[e3] = (lossless 
recording device). For the propagating part of the spectrum [k y < uj/cRe[e^\) however, 
T / for any value of the recording device's permittivity 63, confirming the point made 
in the introduction section that presence of the recording device does not alter the object's 
information content transported by the propagating part of the photon spectrum. Therefore, 
the dissipation in the recording device plays a fundamental role in transferring energy as 
well as fine-structure information about the object (carried by the evanescent components 
of the spectrum). If there is no recording device present, there is no energy or information 
flow in the system for all of the modes with k y > uj/cRe[e^\. As soon as the recording device 
is introduced, the flow of energy and information (for modes with k y > u/cRe[e^\) from the 
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object to IRD starts. Therefore, the dissipation in the recording device in a way materializes 
the fine structure information about an object. 

Moreover, not only does the Eq. [H] determine the amount of energy transferred, but it also 
establishes the amount of information "incident" on the recording device. The total energy 
absorbed in the device can be found by integrating expression [S] over the whole detector's 
length in y direction to give, 

JL = \J \<ky)\ 2 Re\UK)] \B (k y )\ 2 dk y , (10) 

where E/L is the total energy absorbed per unit detector's length and t is the record- 
ing time. It would be useful to define the energy spectral density function S x = 
\ \ T (k y )\ 2 Re[^3(k y )} |-Bo(^y) | 2 , which may be further used to define the density probability 
function for occurrence of k y spatial frequency at the imaging plane, 

P™9{k y ) = M/| |r| 2 ite[£g] \B (k y )\ 2 = *-S x (hy), (11) 

where \I/ is the normalization constant. Using the above equation as well as the definition of 
information (jSJ), one can characterize the amount of object's information I 1 ™ 9 at the imaging 
plane, 

c 9 = -&r 9 io g2 [pn (12) 

Ky 

where pf" 9 = |^ f^gj^ \ T \ 2 ^[£3] \Bo(k y )\ 2 dk y . Next, one may define the function 
Alg = (ig — J* mfl ) (let us call it information distortion function) that quantifies the dif- 
ference between the "true" information content of an object and the transferred information 
content. Clearly, the goal of the imaging system would be to provide maximum energy 
flow between the object and imaging planes, minimizing the information distortion function 
A/^. In reality however, the "true" information content of an object 1$ is unknown to the 
observer (otherwise, what is the point of imaging the object) and minimization procedure 
is not possible. In this situation the only viable approach is to find a set of parameters 
for which the optical transfer function T of the imaging system stays constant over wider 
range of the wave numbers k y . Nonetheless, it would be instructive to see how variation of 
a single variable (dissipation in the recording device) in the imaging system influences both, 
the total energy absorbed in the detector and the information distortion function A/5 for 
the case when the initial object is a two-slit system. Figures 0] and [5] show both quantities 
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versus the imaginary part of the recording device's permittivity 7 (e^ = .Re [€3] +2*7). As 
one can see, there is a non-monotonous dependence of both quantities on the amount of 
dissipation in the recording device. Furthermore, there also seems to be a range of values 
for the dissipation constant where the energy flow and information distortion function reach 
their corresponding maximum and minimum values. This indicates the existence of some 
advantageous device operation regimes, warranting future investigations. In addition, figure 
6 shows the x component of time-averaged Poynting vector as a function of y coordinate, 
(calculated at the position of the recording device) for the case of two slits imaged by light 
with wavelength 0.52 fim. As one can see, the actual image "imprinted" in the device is 
somewhat different from that shown in Figure [2J, that was calculated using only the transfer 
function for the magnetic field. 

Using the universal bound on the energy cost per bit of information transferred, one 
arrives at the following inequality for the recording time, 

f 2 > ^^pjims (13) 
TlL 

Since both, the normalization constant \I/ and the amount of information available at the 
imaging plane P™ 9 are functions of the recording device's dissipation constant, the above 
inequality can in general be written as the relation between the recording time and the 
dissipation constant of the recording device (the inequality would not be universal but 
object specific). 



V. SUMMARY 



In conclusion, we have shown that the recording device plays a fundamental role in the 
image formation and registration. Not only does it change the optical transfer function of 
the imaging system, but it actually determines what the observer will see of the imaged 
object, depending on the physical characteristics of the recorder itself. In the case studied 
here, the information contained in the observed object depends on the imaging system's 
characteristics, and the question of how much the recording device affects the result of the 
observation becomes important. The critical influence of the recording device comes from 
the requirement of the presence of both (incident and reflected) evanescent waves to provide 
a finite energy flux between the object and the recording device. As a result, the issue of the 
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properties of the recording device becomes non-trivial for the superlensing where a significant 
information is transmitted in the evanescent part of the spectrum. On a more philosophical 
note, we could not help noticing that the results of this study have certain parallels with the 
"measurement problem" in quantum mechanics where the observed property of an object 
emerges from the actual act of measurement. 
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FIG. 1: A schematic geometry of a system under investigation. 
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FIG. 2: The image of two slits of width 2.5 nm, separated by 20 nm. The solid line represents the 
case when there is no recording device present (Ime3=0) and the dashed line represents the situation 
with the detector present (e3=10. 47+8. 0*1). The incident light wavelength is 0.52 fim, €q = 11.7, 
ei=-10.47+0. 2811*1, e 2 = 10.47. The thickness of the lens is k d = 0.07 where fc =1.2079x 10 7 
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FIG. 4: The energy absorbed in the detector (per unit length and time) versus its dissipation 
constant (7 = Ime^) for the case of two slits imaged using light with 0.52 jum wavelength. The 
parameters of the imaging system are the same as those given in Fig. 2 
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FIG. 5: The information distortion function AIg=(Is — I 1 ™ 9 ) versus the dissipation constant 7 of 
the recording device for the case when two slits are imaged using light with 0.52 fim wavelength. 
The parameters of the imaging system are the same as those given in Fig. 2. The degree of 
discretization 5 = 0.01 * ko, where ko = 1.2079 x 10 7 mT 1 . 
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FIG. 6: The calculated time-averaged energy flux (the x component of the Poynting vector) of 
an image of two slits of width 2.5 nm, separated by the distance 20 nm. The parameters of the 
imaging system are the same as those given in Fig. 2 
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